On 
On 

IX, 



> 



QCD 



D. Espriu* 

Departament d'Estructura i Constituents de la Materia and IFAE, 

Universitat de Barcelona, 

Diagonal 647, E-08028 Barcelona. 



Abstract 

This is the written version of the lecture on deep inelastic scattering and related top- 
ics in QCD, delivered in the course of the XXVI International Meeting on Fundamental 
Physics to an audience of young experimentalists. The aim is fundamentally pedagogical. 

£\j ■ I review the theoretical setting of the Altarelli-Parisi equations, discuss recent determi- 

nations of a s from deep-inelastic scattering and then move to the kinematical region 

Q\ • explored by HERA. In the way I mention some unsolved theoretical problems. I discuss 

low- a; physics and to what extent log - resummations are called for. 
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1 Introduction 

These notes are not an introduction to Quantum Chromodynamics (QCD), the theory of 
strong interactions. Many excellent textbooks exist where the interested reader can find 
clear expositions of the subject Q. In fact we shall assume here that the reader is familiar 
with the basic technical tenets of perturbative QCD, such as Feynman diagrams, dimensional 
regularization, renormalization, beta function and so on. We pretend rather to give physical 
insight of the reasons behind a rather remarkable theoretical development which is nearly as 
old as QCD itself, namely deep inelastic scattering. Why is it so remarkable? 

QCD is, in a way, a rather simple theory (specially when compared to the intricacies of the 
electroweak part of the Standard Model). It is just a simple extension of good old Quantum 
Electrodynamics. Instead of matter fields carrying electrical charge +1 (and anticharge -1), 
if we are talking about electrons, the matter fields of QCD (the quarks) carry a new quantum 
number: color. Color can take three different values (and their corresponding anti- values) . 
Furthermore the intermediate bosons (the gluons), unlike photons which cannot change the 
charge of a particle, change the color of a quark. They may for instance turn a red quark 
into a blue quark. This simple fact implies that the gauge group of QCD is much larger that 
the U{\) of QED. Since every quark comes in three copies, whose labels get exchanged, there 
is a SU(3) invar ianceQ. 

Simple as this theory may seem, it is not an easy matter in QCD to relate theory and 
experiment. It is well known that the fields and particles we know how to compute with (with 
the simplest tool at our disposal, perturbation theory) are not those that are observed by 
experimentalists in their detectors due to the phenomenon of confinement. Quarks and gluons 
are real, but they cannot be detected as free particles as they are known to be confined inside 
hadrons. In view of this is quite remarkable that there are theoretical techniques enabling us 
to put the theory to very stringent tests. 

The phenomenon of confinement comes about because the coupling constant of QCD 

(which is relatively small at large values of the momentum transfer, the value quoted by the 

Particle Data GroupQ is a s (Mz) = 0.119 ± 0.002; notice the truly amazing precision, which 

may soon be reduced to a mere 1%) becomes strong as the energy decreases. The behaviour 

1 The reason why the symmetry group is 5(7(3) and not (7(3) — which would have nine gluons — has to 
do with the mathematical fact that e a p 1 is not an invariant tensor of (7(3). This tensor is necessary to build 
combinations of three quarks which are antisymmetric, as required by Fermi statistics. For instance, A ++ = 
\v}u'u'). Being the lightest hadron with this quark contents we expect to have the three quarks in the ground 
state, hence in a symmetric wave function. This is in contradiction with Fermi statistics. The contradiction 
can be solved if we admit the existence of a new quantum number a and \n}u^u^) = -y-e a ^ J \ul,uLuJ f ). 



predicted by the renormalization-group is, at one loop, 



"•«> = fio g (Q7AW (1) 

where Aqcd is a renormalization-group invariant, but scheme dependent, quantity and /3i = 
-11/2 + JV//3. The preferred value is Aqcd = 219±|| MeV (5 flavours, MS-scheme). The 
meaning of the scale-dependent, renormalized coupling constant is roughly the following: it 
is the "effective" coupling, relevant at the scale Q, namely, the one that (within the choosen 
scheme) minimizes further quantum corrections, in particular resumming all large logs. From 
([[]) we see that precisely at the scale Q 2 = Aq CD , the effective coupling has a pole. Of course 
well before that scale is reached perturbation theory becomes completely unreliable, and the 
1/r potential of the perturbative interaction change to a stronger behaviour, possibly to a 
linear ~ r behaviour. When computing the production of any physical hadron (typically of 
mass ~ few Aqcd), perturbation theory is completely useless. 

One instance where perturbative QCD can be applied is to inclusive processes, provided 
that the characteristic momentum transfer is large enough. These will not be discussed in 
this lecture. The interested reader can look at the determination of a s through Rhad Rt, f° r 
instance in.[[| 

A clear application of perturbative QCD is provided by deep inelastic scattering. The 
subject is now over twenty years old and, by now, perturbative QCD has been tested to a 
high degree. Furthermore, the commissioning of HERA has opened a new kinematical region 
where it may be possible to study the onset of non-perturbative effects in a controlled fashion. 
The exploration of this region is a fascinating subject interesting on its own right. 

Due to the lack of time and space we have not included two sections that, in our view, 
should be in any general review of perturbative QCD and deep inelastic scattering. The 
first one concerns the so-called "spin of the proton" problem Q]. Another topic that is not 
covered at all is the study of the photon structure functions. M The list of references is very 
incomplete and those provided merely reflect personal tastes. 

2 Logs in QCD 

Beyond tree level most Feynman diagrams are ultraviolet divergent. Take for instance the 
one diagram contributing at one loop to the gluon propagator. Neglecting external momenta, 
the integral over the momenta of the internal particles is of the form 

d 4 k k a k? 



(2tt) 4 k 4 



oo. (2) 



To make sense of the theory and get a finite result we must introduce a cut-off A and 
counterterms. A possible method is to perform a subtraction at some q 2 = — /i 2 . For instance, 
for the self-energy of the gluon propagator 

U(q 2 ) - U(-fi 2 ) = U R (q 2 ) = finite. (3) 

Alternatively we can make sense of the integrals using dimensional regularization by 
continuing the dimensionality from 4 to n = 4 + 2e, 

d A k f d n k 



(2vr) 4 ' J (2vr)^' ^ 

and subtract just the poles in 1/e (minimal subtraction, MS) or also the 7^ — log4-7r that 



always accompanies the singularity in 1/e (improved minimal subtraction, MS). For in- 
stance, for the quark contribution to the gluon self-energy, one has the following result after 
computing the integral in n = 4 + 2e dimensions 

n(<? 2 ) = -^5 ab {\ + 1E + log ^_ + ...). (5 ) 



Using the MS and MS schemes one gets 



1 1 us(<r > - --^Ui ~t + log -^ + . . .) iWg 2 ) = -^Mlog K + ■ ■ •) ( 6 ) 



■^(7 B + log ^ + ■ • n^ (g ) = -^^(bg -^ 



The above expressions illustrate the appeareance of ultraviolet logs logq 2 /^ 2 through the 
renormalization procedure. There are, however, other source of logs in QCD. They are of the 
form logc/ 2 /A 2 where A 2 can is some external momentum squared or a small (mass) 2 that we 
have given by hand to the (massless) gluon. While the former are associated to ultraviolet 
divergent integrals (integrals with a bad behaviour when the internal momentum is large), 
the latter are infrared logs and are related to Feynman diagrams with a bad behaviour when 
one or more external momenta vanish. Ultraviolet logs appear in any renormalizable field 
theory after renormalization. On the contrary, infrared logs appear whenever a theory has 
massless particles in the spectrum (such as photons or gluons). A given Feyman diagram can 
give rise to both type of singularities at the same time. 

There actually two classes of infrared logs caused by massless particles. The so-called 
infrared divergences arise from the presence of a soft massless particle (k^ — > 0). For instance 
in the process e + e~ — > /x + /i~ at the one loop level we have to compute the integral (figure [j]) 

(7) 



(27r) 4 k 2 [(pi + k) 2 - m 2 ][(p 2 + k) 2 -m 
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Figure 1: Example of diagram leading to an infrared divergence. 



When p\ = p\ = rn 2 the integral behaves for k^ — ► as 



d 4 k 1 



(8) 



(2vr) 4 A: 4 ' 

and diverges. This divergence is unphysical so it must be cancelled by something else. The 
Bloch-Nordsieck theorem || states that in inclusive enough cross-sections the infrared logs 
cancel. What do we mean by 'inclusive enough'? A detector will not be able to discern a 
'true' muon from a muon accompanied by a soft enough photon (with k — ► 0). Therefore, we 
have to consider diagrams where a soft photon is radiated by the muon, square the modulus 
of the amplitude and integrate over the available phase space (which actually depends on the 
experimental cut). When this is done the result is infrared finite. The relevant diagrams are 
depicted in figure g. 





Figure 2: Real and virtual photons have to be included for IR safe results. 

The other type of infrared logs are called mass singularities. They occur in theories with 
massless particles because two parallel massless particles have an invariant mass equal to zero 



k 2 = {k 1 + k 2 f = ||(wi+«2,0,0,wi +UJ2) 



0. 



(9) 



The appeareance of such a mass singularity is illustrated in figure 3 

q 



Figure 3: Diagram with a mass singularity. 



1 _ 1 

(p - k) 2 ~ p 2 + k 2 - 2k°p° + 2k°p° cos 9 ' ^ 

the denominator vanishes when we set all particles on shell (p 2 = k 2 = 0) and 6 — > (i.e. fe 
is parallel to p). Even if one of the two particles is massive there is a singularity, provided 
the 3-momenta are parallel. 

The Kinoshita-Lee-Nauenberg theoremj?]] ensures that for inclusive enough cross section 
the mass singularities also cancel. Both for mass singularities and for infrared divergences 
there is a trade-off between A 2 , the infrared regulator of a massless particle, and the energy 
and angle resolution of the inclusive cross section AE, A9. 

In practice, it is better to regulate the infrared logs using dimensional regularization 
(introducing A 2 leads to difficulties with gauge invariance). Real gluon emission diagrams are 
regulated by performing the phase space integration in n dimensions. 

There is in fact a lot of physical insight hidden in the infrared logs. Physical arguments 
tell us that the probability of finding a 'bare' isolated muon should be zero. We know this 
because detectors are unable to tell apart a muon from a muon plus one soft photon or indeed 
from a muon plus any number of soft photons. Infrared divergences in QED can be summed 
up and then one sees that the probability of finding an isolated muon is indeed zero and not 
infinite as the one loop diagram led us to believe. Whenever a Feynman diagram is infrared 
divergent it means that we have forgotten something relevant. 

Let us consider in QED the interaction of a charged fermion with an external source and 
let us expand in the number of virtual photons n (i.e. in the number of loops). The total 
amplitude will be expressed as 



M(p,p') = J2M n (p,p') (11) 



n=0 



then a calculation shows that 



M = m , 

Mi = rriQaB + m\, 

Mi = mo \- miaB + rri2, 



(12) 



The quantities m n are IR-finite, while B is IR-divergent. The series in ( |Tl"D can be summed 
up 

oo 

M = ex.p(aB) J^ m n , m n ~ a n , (13) 

n=0 
6 



and B can be obtained just from the lowest order diagram. Introducing an IR cut-off A, 
B ~ — logm 2 /A 2 , which indeed shows that when we remove the cut-off the probability of 
finding an isolated charged fermion is zero in QED. The addition of soft photons changes that 
result, multiplying the total amplitude by a factor ~ (AE/X) 2 . There is a trade between the 
infrared regulator and AE,A6. The latter are, of course, detector-dependent quantities. 

Although only partial results exist ||, it is believed that a similar exponentiation takes 
place in QCD. Due to the confinement subtleties it is unclear whether the suppression factor 
is compensated by radiation of soft gluons. Even if this compensation does actually take 
place that would not disprove confinement, only that confinement would have nothing to do 
with the structure of infrared singularities of the theory. 

The previous discussion can be summarized in the following way: due to IR singularities 
one is forced to consider cross sections not of individual particles in the final state, but rather 
of bunches of particles, each 'hard' quark and gluon surrounded by a 'soft' cloud of gluons 
and, perhaps, quarks. These bunches are called 'jets'. 

The Bloch-Nordsieck and Kinoshita-Lee-Nauenberg theorems guarantee the finiteness of 
the cross-sections. We have to define an energy and angle resolution. For instance, if p is 
the momentum of a primary quark we can impose that the energy of each soft particle in 
its jet satisfies k® < epo and also that arg(p, ki) < 5. We will get singularities of the form 
a s log e log 8 when e, 5 — ► 0. The specific details depend on the precise definition of the jet. 

3 Free Parton Model 

The counterpart of having an effective coupling constant which grows at low energies is that 
the theory becomes simple at high energies, making perturbative calculations possible (at least 
some of them). Indeed a brilliant confirmation of the existence of nearly free constituents 
inside the nucleon was provided more than twenty years ago by a series of experiments 
carried out at SLACpl. Then it became possible to scatter electrons off nucleons in fixed 
target experiments with a typical momentum transfer ~ 1 — 10 (GeV) 2 , a kinematical range 
unexplored until that time. The kinematics of Deep Inelastic Scattering (DIS) processes is 
shown in fig. 4. 

The virtual intermediate boson is far off its mass-shell and scatters off a quark or gluon 



in a time of 0(1/ y/— q 2 ). Typically quarks and gluons are themselves off-shell by an amount 
of O(Aqcd)- After the scattering the outgoing particles recombine into hadrons in a time of 
0(1/Aqcd)- Thus Deep Inelastic is a two-step process 




Figure 4: The kinematics of deep inelastic scattering. 

• Short distance scattering occurs with a large momentum transfer. Well described by 
perturbation theory. 

• Outgoing particles recombine. Not calculable in perturbation theory. 

However, the second part can be side-stepped all together for fully inclusive rates. Then 
perturbation theory is adecuate to describe many features of DIS. 

If we place ourselves in the center of mass of the hadron and virtual intermediate boson 
both particles move very fast towards each other. Whatever components the hadron contains 
they will all have moments parallel to P^, up to transversal motion of O(Aqcd)- Let us 
write 

p v = x pf. (14) 

The squared CM energy of the lepton and proton constituent will be 

s = ( x p + k) 2 ~ 2xPk ~ xs. (15) 

We neglect masses (as well as the fact that constituents are off-shell by O(Aqcd)) The final 
momentum of the constituent is xP + q. Therefore 

0~ (xP + q) 2 ~2xPq + q 2 , (16) 

so x = —q 2 /2Pq. If v is the energy transfer in the LAB system, we can also write 

-q 2 
x = — ^— (17) 

zuniN 

771,7V being the nucleon mass. It is convenient to introduce 

Pq Phf , s 

In the lab frame y = v/E and < y < 1. y is thus the relative energy loss of the colliding 
lepton. 



Let us for the time being ignore altogether QCD interactions and let us assume that 
constituents of the nucleons (which we will call partons) are free. DIS will then be described 
by an incoherent sum over elementary processes. The partonic differential cross sections in 
the LAB frame will be 



vq, i/g-scattering 



~fi = ^ ) ( g2 -M^ + SnO- ~V)]> ( 19 ) 

dot g 2 2 irmE 2 2 2 



eg-scattering 



da e 2 ^a 2 mE 

dy q 4 



^ = Q 2 ^[l + (l-y) 2 ]. (21) 



The neutral current sector is dominated by 7 exchange below q 2 = M§, so we have not 
bothered to include Z exchange. In (^) Q is the quark electric charge (in units of e) and m 
is the target mass. Since p^ = xP^, we just take m = xm^. Then, for instance, 

d 2 a e 2 \-ko?: 

dxdy " q A 



Q 2 ^^[l + (l-y) 2 ]. (22) 



Let u(x)dx, d(x)dx, ... be the number of u, d, ... quarks with momentum fraction between 
x and x + dx in a nucleon. Then xu(x),xd(x), ... will be the fraction of the nucleon momentum 
carried by u,d, ... quarks. We, of course, identify quarks with partons and, since we assume 
that they are free, proceed to sum incoherently over the different scattering possibilities. For 
instance in ep — > eX 

d 2 a 2ira 2 1 + (1 — y) 2 A . , . ... 1 . „ . -,. .. 1. . . , .., . . 

— - = L_«y_ ( U x ) +u(x)) + -(d(x) + d (x ) + - (s (x ) + s(x ))]. (23) 

dxdy s xy L 9 9 9 

(We neglect here the possible contribution from the sea of heavy quarks in the nucleon.) 
Other DIS processes weigh differently quarks and antiquarks. For instance, in up — * fiX if 
f 2 

d 2 a G 2 F s 



-q 2 -C Ml we have 



x— b —[c 2 c d(x) + s 2 c s(x)+u(x)(l - y)\ (24) 



dxdy ir 

with c c = cos# c , s c = sin^ c , the cosinus and sinus of the Cabibbo angle, respectively. 

The parton distribution functions (PDF) q{x) are quantities which are not calculable 
within perturbative QCD, as we will see. 



Probably the first thing that one learns is that gluons are very important. From the 
SLAC-MIT data§ 

Q = U + D + S = J dxx(u(x) + d{x) + s(x)) ~ 0.44, (25) 

Jo 

Q = U + D + S = [ dxx(u(x) + d{x) + s{x)) ~ 0.07. (26) 

Jo 

The total fraction of momentum carried by quarks (and antiquarks) is only about 50% The 
rest is carried by gluons (parametrized by a PDF g(x)), showing that although the naive 
quark model works very well is just a gross simplification as a model of hadrons, at least at 
large — q 2 . In fact we know the asymptotic values of ( |25| ) and ( |26| ) based in the equipartition 
of energy in a free theory (since, asymptotically, QCD is free) 

dxxq(x) — > — — / dxxq(x) — ► — -. (27) 

o v J 16 + 37V/ Jo yy J 16 + 3N f v ; 

From the above limiting values we see and at higher energies the total momentum carried 
by constituent or valence quarks diminishes and that an equally important role is played by 
particles from the Dirac sea of the nucleon. 

Another example where the quark model fails to describe some basic features of hadrons 
is provided by the 'spin of the proton' problem Q]. /^-scattering on polarized targets shows 
that the fraction of the total spin of the proton that can naively be associated to constituent 
quarks is surprisingly small. We shall not dwelve on this matter further here. 

Nevertheless, there are some obvious sum rules for the parton distribution functions which 

can ultimately be explained in terms of the quark model. For the proton 

-I 

dx(u(x) - u(x)) = 2, (28) 

o 

l 

dx(d(x) - d{x)) = 1, (29) 

o 

i 

dx(s(x) - s{x)) = 0. (30) 

o 

On QCD grounds we expect that this free parton model description of the hadrons becomes 
more and more accurate when — q 2 — > oo, v — > oo, while keeping x fixed. This limit is known 
as Bjorken scaling and in the strict —q 2 = oo limit everything depends just on x. 

Let us now try to rederive the previous results in a more theoretical setting. Let us 
consider for instance vp scattering. Then 

d 2 a _ G 2 F m N jiiVjJ 
d{—q z )dv ns A 
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where 

IT = ±Tr[f(l - TshVU " Kh^kakp (32) 

is the trace over the leptonic external lines, and H^ u is given by 

J2(P\M0)\X(P'))(X(P')\J U (0)\P) = fd 4 ze^(P\J,(z)J u (0)\P), (33) 

x J 

which is just Imn^;,^), with 

U^(q) = |d 4 z e ^(P|TJ^z)J,(0)|P). (34) 

We decompose H^ u as 

H^ = -g^Fi + ^F 2 + -5- W P'/F 3 (35) 

(If we assume that we are working with non-polarized targets P and q are the only vectors at 
our disposal.) i*\, F 2 and -F3 are called the nucleon structure functions. Using the kinematical 
relations x = — q 2 /2vniN and y = 2vn^vjs we get 

d{— q )dv = usdxdy, (36) 

(]1 ° ^ [iW + F 2 (l - y) - F 3 xy(l - V -)\. (37) 



dxdy 27r " " 2' 

Let us now compare with the free parton model. We see that (restoring the up index, to 
make apparent that the structure functions are process dependent) 

Ff (x) = c 2 c (u(x) + d{x)) + s 2 c (s(x) + u(x)), (38) 

F 2 up (x) = 2xc 2 c (u(x) + d(x)) + 2xs 2 c (s(x) + u(x)), (39) 

F% p (x) = 2c 2 c {u(x) - d{x)) + 2s 2 c {-s{x) + u{x)). (40) 

For other processes the actual expressions may vary but the structure functions are always 
linear combinations of the parton distribution functions, i.e. F 2 (x) = xJ2fdfQf( x )i e ^ c - 
Note that in the free parton model 

F L {x) = F 2 (x) - ^M = o. (41) 

This is the Callan-Gross relation, which actually is not an exact one; it gets modified when 
the q 2 dependence is included, i.e. we depart from the strict — q 2 = oo limit. 
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An exact sum rule, which is easily expressed in terms of the structure function i^fa) was 
given by Adler 

/ -(Ft 1 -Ff)=4(I 3 ) (42) 

where ^3 is the third component of the target of isospin /. Other sum rules are not exact, 
except in the strict free parton model, but their violations are computable within perturbative 
QCD. Two instances are the Gross-Llewellyn-Smith sum rule 

- / dx(F%> + F* p ) = f dx{u(x) - u(x) + d(x) - d{x)) + 0(a s ) = 3 + 0(a s ), (43) 
2 jo Jo 

and the Gottfried sum rule 

f 1 dr 1 f 1 _ 

/ — (FF - FD = ^ + 0(a s ) + / dx(u(x) - d(x)). (44) 

Jo x 3 Jo 

Violations to the different sum rules are a theoretically clean way to extract a s . Of course 
in practice things are difficult because the sum rules involve an integral over all values of x, 
which is always poorly known in some range of the integrand and extrapolations are needed. 
If we assume that a s is extracted from some other source, the Gottfried sum rule provides 
some interesting information on the sea contents of light antiquaries in the nucleon. The 



collaboration NMC[10| has determined that for the proton 

U-D= [ dx(u(x) - d{x)) = -0.15 ± 0.036. (45) 

Jo 

In the proton quark sea there are many more d-type antiquarks than it-type. This is in fact 
a large isospin violation, much larger than expected to the mass difference of the quarks and, 
in fact, goes, at least naively, in the opposite direction. 

The isospin violation is larger at low values of x {x < 0.2) and also at low values of 
Q 2 , which hints that long-distance physics (quite remote from perturbative QCD) is called 
for. The enhancement of d-type antiquarks is confirmed by other experiments. For instance 
Na5l|ll|] finds at x = 0.18 u/d = 0.51 ± 0.09, while NuSeaC||] finds at Q 2 = 7.4 GeV 2 that 



U -D = -0.10 ±0.024. 

Could this be due to the exclusion principle, as illustrated in figure 5, which makes harder 
for u-type antiquarks to appear in the sea? It could be, but it is very difficult to come up with 
quantitative results. A partial understanding is provided by the chiral quark model, including 
pion exchange (see figure 5), but this type of physics is still very poorly understood. 
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Figure 5: Valence and sea quarks. Contribution from the chiral quark model. 

4 Scaling Violations 

It is plain clear from the data that there is some Q 2 = —q 2 dependence in the structure 
functions. In other words, there are violations of Bjorken scaling and actually Fj = Ff(x, Q 2 ). 
The free parton model is not completely correct (no big surprise, of course). Our job is to 
try to understand these violations in the framework of QCD. 

Let us assume that we have isolated a parton with initial momentum p = xP. The 
probability of finding such a parton is given by q(x). At the parton level the structure 
function F<i is just F2 = x5f (<5/ is the appropriate charge. At the proton level, however, 
this partonic cross-section has to be multiplied by the probability of finding the parton with 
momentum fraction x, i.e. by q(x). We write this in the form 

F 2 (x) = xS f I' d$q(Q6(x -0 = xS f C ^?(0*(1 " ?)■ (46) 

Jo Jo 4 4 






Figure 6: Contribution at 0(a s ) to the relevant DIS subprocess. 

At 0{a s ) many diagrams contribute. They are given in figure 6. We are looking for 
scaling violations and, therefore, we must look for logs. In other words we must investigate 
ultraviolet, infrared and mass singularities of any kind. It turns out that ultraviolet singu- 
larities are proportional to the free result and simply renormalize 5f. Infrared divergences 
cancel amongst all diagrams. Only the mass singularity present in diagram (d) when the 
momentum of the gluon is parallel to that of the gluon survives. Keeping only logarithmic 
terms, the caculation at 0{a s ) amounts to the replacement 

*(l-f)-*(l-f) + gj°(f)log^ (47) 
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where A is an infrared regulator and 



I + z 2 3 



P(z) = C p\.wt^r + 2 6{1 ~ z)] - (48) 

Only the logarithmic term is retained for this discussion. Now we understand the reason for 
writing things in apparently such a complicated way. First of all, scaling violations appear 
through the contribution of real soft particles. £ is the original momentum fraction of the 
nucleon carried by the parton, which is reduced to x < £ after the emission of the soft gluon. 
The hard scattering takes place with the parton carrying fraction x. All along the discussion, 
the transverse motion of the partons inside the target is neglected as well as are all masses. 
Then 

F 2 {x)=x8 f [q (x)+£^q (t)£p(Z)] oe &]. (49) 

In addition we have replaced q(x) by qo(x), the bare PDF. If we now define the renormalized 
PDF by 



j(x, li l ) = q (x) + ^J x jqo(x)P(^) log jj, (50) 



we can write 

F 2 (x,Q 2 )=x6 f [q(x, fJ 2 ) + f l ^q( x ,^p£)\og%}. (51) 

No doubt the similarity with the usual renormalization process did not go unnoticed. Now the 
infrared regulator has been eliminated (hidden in the bare PDF) at the expense of introducing 
a renormalization-scale dependence. These are the sought after scaling violations. 

5 Altarelli-Parisi Equations and Aqcd 

At this point it is convenient to introduce the variable t = h log/i 2 /Ag C - D . It then follows 

from © that 

d , . a s (t) f 1 df .. XT ^,x, ,„„. 

-q(x,t) = -^- / -j-q (tt )P - , 52 

at it Jx £ i 

which immediately translate into differential equations for the structure functions themselves. 
These are the Altarelli-Parisi equations^]. They summarize the rate of change of the parton 
distribution functions with t. 

We define the moments of the PDF's by 

q(n,t)= [ dxx n ~ l q{x,t). (53) 



/o 

Introducing the anomalous dimension 7„ as 



7n == / dxx"~ 1 P(x), (54) 

Jo 
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the convolution over the fractional momentum £ transforms into a product 

d &s(t) 

-£;q(n,t) = -y n q(n,t). 

at tt 

This leads to following scaling behaviour for the moments of the structure functions 

— 

F 2 (n, Q ) = F 2 (n, Q ) — — 

\ a s (Q) 



(55) 



(56) 



which is our final expression. Experiments agree on the whole very nicely with the scaling 
violations predicted by QCD. Taking into account all the subtle points of Quantum Field 
Theory that have gone into the analysis, this provides a beautiful check of the theoretical 
framework. 

We have been considering F 2 , but the same procedure can be repeated for any structure 
function. The expression ([56]) amounts to resumming the leading logs obtained by iteration 
of soft collinear gluons. The diagram is the one shown in figure 7. 




Figure 7: Handbag and ladder diagrams. 

As a simplifying hypothesis we have neglected mixing. In fact, the evolution equation is 
a (2Nf + 1) x (2Nf + 1) matrix, involving quarks and gluons. In the flavour singlet case life is 
more complicated; there is mixing with gluon operators and therefore one must also consider 
gluon parton distribution functions as well 

dq(x,t) a s (t) f 1 dy . 



Of 

dg{x,t) 
dt 



IT 



.q(y,t)P qg (-) + g(y,t)P gq (-)] 

x y y y 



^ r%{y,t)P gg {-) + q{y,t)P qg £) 

tt Jx y y y 



(57) 



(58) 



The detailed form of the Altarelli-Parisi kernels at leading and NLO order can be found in 
[14|. No complete calculation exists yet at the NNLO to my knowledge, just some partial 
results. 

It is important to realize that the Altarelli-Parisi equations are not exact. They take 
into account the perturbative contribution only (and this up to a given order in perturbation 
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theory). They also neglect transverse motion, which leads to corrections of 0(Aqcd/Q 2 ) 
to the leading results. These are more easily dealt with in the perhaps more rigorous (but 
more cumbersome) treatment based in the Operator Product Expansion fi~5|| , which will not 
be discussed here. Target mass corrections should be equally taken into account. They 
are particularly important near thresholds (such as the charm and bottom thresholds). For 
instance the proper treatment of thresholds is highly relevant for HERA, since a big chunk 
of the data comes from a region close to these thresholds. And, of course, mass corrections 
are important for charm PDF from the sea, which have actually been recently measured. 

The analysis of Deep Inelastic Scattering based on the Altarelli-Parisi equations (or, 
alternatively, on the Operator Product Expansion) has been one of the most clear tests of 
perturbative QCD and traditionally the best way of determining a s , which, as we have seen, 
enters in the scaling violations. However, at present the value of a s {Mz) extracted from 
Z-physics is equally accurate if not more. 

For some time a discrepancy was claimed between the value of a s obtained from the 
analysis of scaling violations and the Z-pole value. For instance the value quoted by CCFR 
was a s (Mz) = 0.111 ± 0.004, well below the world average[||, and quite away, for instance, 
from measurements based on event shapes at the Z peak (a s (Mz) = 0.122 ± 0.007; for 



measurements of a s at LEP, see e.g. [|17|1 ) ■ The average a s from DIS was given[18] just two 
years ago to be a s (Mz) = 0.113 ± 0.005. Theoretical speculations were fuelled and it was 
claimed that non-perturbative corrections could be larger than originally thought. 

Deep inelastic scattering data have been reanalyzed recently and the quoted value for 



a s (M z ) from DIS from a global fit[|T9| to all data is a s (M z ) = 0.118 ± 0.005. The agreement 
with other determinations is now almost perfect. The discrepancy was apparently due, it is 
claimed, to the energy calibration of the detector (in the case of CCFR, at least; the new 
CCFR value is a s (Mz) = 0.119±0.005[|l6|]), a better understanding of higher twist corrections 
(not discussed here) and from the treatment of the contribution from heavy quarks from the 
sea, charm in particular. We will see later, when we discuss in somewhat more detail the 
form of the PDF, that one must actually make a number of hypothesis before being able to 
extract a s (Mz)- While do not regard the issue as totally settled yet, because some of the 
older data sets are very poorly described by the now preferred value of a s (Mz), the most 
recent data coming from HERAJ20J (HI and ZEUS) give values which agree nicely with each 
other and in fact fall well in the high a s range, almost on top of the new average value. These 
ongoing experiments will become statistically more and more significant in the near future 
in the determination of a s (Mz) as more and more data points pile up. 
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6 Parton Distribution Functions 

We do not know in general how to compute the parton distribution functions, even for 
— q 2 — ► oo. Only their evolution can be reliably computed either through the Operator 
Product Expansion of the use of the Altarelli-Parisi equations and this for large enough 
values of — q 2 . The scaling behaviour is governed by the anomalous dimensions. At leading 
order they are 

7„0") = **["£ + -^ + (j + ^ J + 2) - g i] - 2 -^T R , (62) 

where Cf,Tr and (7,4 are group-theoretical factors. 

An interesting issue is the behaviour of the parton distribution functions at the endpoints 
x = and x = 1. The large n behaviour of the moments probes the x — > 1 region. Since 
it is natural to expect that at the kinematical boundaries the parton distribution functions 
vanish, one can make the following ansatz for x — > 1 

q(x, Q 2 ) ~ A(Q 2 )(l - x )KMQ 2 ))-i. ( 63 ) 

Demanding that eq. (|63| ) fulfills the g 2 evolution equation leads to 

r(l + i/(a s (C^))) 33 - 2iVf 



1 f\ Q 

d0= 33-^V> ( 4- 7s) - (65) 



Likewise, for the gluons we have 



^,Q)-^o r(2 + l/(as(Q2))) lQg(1 _ x) • (66) 

The constants Ao,Aq and vq are not calculable on perturbative QCD and depend on the 
specific operator, do is universal. 

When x — ► 1 the gluon distribution functions approach zero more rapidly than the quark 
ones. For large values of x the quark contents of nucleons is the relevant one. Second order 
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corrections to this asymptotic behaviour can be derived in a similar way and are known. It 
turns out that the correction is arbitraryly large if one gets sufficiently close to x = 1. This is 
because the collinear gluon is, in addition, soft in that exceptional configuration thus giving 
rise to a log(l — x) singularity Multiple emission is then kinematically favoured, since the 
log overcomes the a s suppression. 

For small values of x the opposite behaviour takes place, the gluon distribution function 
eventually becomes dominant. At LHC the cross-section will be greatly dominated by low- 
x physics and the important process there will be gluon-gluon scattering. At the current 
Tevatron run the quark contents of protons and antiprotons is still dominant. Let us see why 
gluons dominate completely at low x. 

The key point is the appeareance of the singularity for j = 1. Indeed, as j — > 1 

2iV 
lM ~ — T- ( 67 ) 

Then 

9(ji *) = 9(j, to) exp[— — — log — ]. (68) 

TTPiU — i; to 

Then we proceed to evaluate g(x, t) by performing an inverse Mellin transform 

xg(x,t) = — j djx l ~ j g{j,t). (69) 

2-iri Jc 

The integration circuit is a line in the direction of the imaginary axis in the complex j plane, 
to the right of the j = 1 singularity. The saddle point method can now be used provided that 
log(l/x) is large and this is the reason why this procedure gives only the small x behaviour. 
Working things out we see that the gluon parton distribution function for low x behaves as 



g(x) r, - exp J C(Q 2 ) log-, (70) 

XV X 

where C{Q 2 ) is calculable. Unfortunately, this answer is not totally satisfactory because 
something must stop the growth in g{x) for low x, or else one runs into unitarity problems 
sooner or later, and thus eq. (|70| ) it is not credible all the way to x = 0. Technically 
speaking, there must be corrections that destabilize the saddle point solution. Physically, 
the uncontrolled growth of the gluon distribution is an infrared unstability. The density of 
soft gluons is too large. Shadowing and non-linear evolution equations are the buzzwords 



here [21 



The double scaling limit (high Q 2 , low x) is well supported by the data. See for instance[| 
for a recent analysis coming from measurements of Fi at HI. 
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Except in these two limiting cases, PDF's have to be parametrized. The way one proceeds 
is by proposing a given parametrization at some reference value Qq, then evolve to all desired 
values of Q 2 using the Altarelli-Parisi equations, then perform a global fit of the parameters 
describing the PDF and, at the same time, determine a s . 

7 Confinement 

Aqcd sets a natural scale in the theory. Well above A.qcd perturbation theory makes sense. 
Of course perturbative QCD at large enough energies describes a world of quasi-free quarks, 
interacting with Coulomb-like forces. We know very well that hadronic physics is a very 
different world where quarks are confined into colorless hadrons. As soon as q 2 ~ ^qcd 
perturbation theory is unreliable. It simply cannot explain confinement. 

What confinement means is that there is a force between quarks that does not decrease 
with distance. There is indeed phenomenological evidence (which is supported by lattice 
analysis) that the interquark potential at large distances in QCD is of the form 

V{r)~ak 2 QCD r- b - + ... (71) 

The first term is a confining quark potential. The constant a has to be ~ 1 because Aq CD 
is the only dimensional quantity at our disposal. The Coulombic part is called the Liischer 
term and plays a crucial role in heavy quark spectroscopy[23|. 



At short distances the behaviour of the interquark potential is totally different. It is 
Coulomb-like 

V(r) ~ -2!. (72) 

r 

Different quarks probe the different regimes. Indeed, because a s (rnt) is so small (say 
~ 0.1), the corresponding Bohr radius is r$ ~ 10~ 2 fm, much smaller than Kq CD . The 
coulombic part of the interquark potential largely dominates. (At such short distances the 
linearly rising potential is not at work, the leading confinement effects are ~ r 3 , as discussed 
by Leutwyler some time ago[p4|], but they can be safely neglected at first approximation.) 

Bottom and charm are in a somewhat intermediate position. a s (77ife) is still relatively 
small. The Bohr radius is 10 -1 fm, smaller but comparable to Kq CD . Spectroscopy is 
basically perturbative, at least for the lowest levels, but some non-perturbative effects are 
visible. Charm is really no-man's land. Both perturbative and non-perturbative effects 
compete even for the ground state n = 1. For light quarks the Bohr radius is several fm and 
the confining potential is fully at work. 
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The existence of a confining potential leads to very large multiplicities and jets. One can 
imagine a quark-antiquark being formed at the primary vertex then moving apart. Part of 
their kinetic energy is deposited in the interquark potential as they move away. Very quickly a 
separation r m is reached where the energy deposited is enough to form a new quark-antiquark 
pair, 

^QCD r m ^ 2m g , (73) 

at that moment the quark-antiquark 'string' breaks and the process is repeated until the 
average relative momentum is small enough and hadronization takes place. 

There is a lot of physics in the string picture. We can think of color forces being confined 
in some sort of tube or string joining the two moving quarks. The chromodynamic energy is 
thus stored in a relatively small region of space-time. If this picture is correct we should expect 
hadronization to take place in this region in preference to any other. This is indeed the case; 
in three jet events (which originate from qqg, with a hard gluon) there is a clear enhancement 
of soft gluon and hadron production in the regions between color lines (representing the 
gluon by a double color line, or qq state), and a relative depletion in other regions. This 
phenomenon is called color coherence J~~ 

8 Dual Models 

We now backtrack in history to the pre-QCD days, and recall that in the 60's the duality 
hypothesis was much in fashion. The hypothesis stated that in strong interactions a sum 
over intermediate states in the s-channel should reproduce the sum over resonances in the 
t-channel. Mathematically, 

Of course for this to have even a chance of being true an infinite number of intermediate 
states is required. It should be stated right away that the evidence for this peculiar property 
was (and still is) rather weak. 

However, in 1968, Veneziano[p6| took the idea seriously and proposed the following am- 
plitude 

. . . T(—a(s))T(—a(t)) , . , , , ,__. 

A ( s ' * = rv \ 1 . «W = « + a's. 75 

1 {—a{s) — a{t)) 

This amplitude is manifestly dual. In 1969 and 1970 Y.Nambu and others unveiled the 

relation between the Veneziano amplitude and open string theory. It was later generalized 

to closed strings by Koba and Nielsen. 
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This of course is the way to make contact with the long-distance properties of QCD 
that we have discussed in the previous section. If in some kinematical regime QCD can be 
described by some type of string theory, an amplitude of the Veneziano type should describe 
strong interactions in a regime where perturbation theory is not valid. 

Unfortunately life is not so easy. First of all, consistency of string theory requires a(0) = 
1 and then the Veneziano amplitude exhibits poles in the s-channel whenever s = (n — 
l)/a', n = 0, 1, .... There is a tachyonic scalar/pseudoscalar particle and a massless vector 
particle. In addition, the amplitude does not exhibit the proper chiral behaviour (Adler 
zero), i.e. ^4(0,0) = 0, if we are to interpret the pseudoscalars as pions. It is a complete 
phenomenological failure. 

Lovelace and Shapiro |27]] finally proposed an amplitude with the correct behaviour or, at 
least, not manifestly incorrect. It is inspired from the theory of supersymmetric strings and 
the equivalent of the Veneziano amplitude now reads 

r(l-q( S ))r(l-q(t)) 
A{8 > t)= T(l - a(s) - a(t)) ' ^s) = a(0) + as, (76) 

where a(0) in principle also equals 1 here. Again this amplitude (which is tachyon-free - 
this is the main virtue of the supersymmetric string) is not physically acceptable because the 
amplitude does not have the correct Adler zero. However, it does reproduce the right chiral 
behaviour if one replaces by hand this value by 1/2 and is then appropriate to describe pion 
scattering. Unfortunately this replacement cannot be derived at present from any known 
string theory. The corresponding trajectory is called the Regge trajectory and corresponds 
to the exchange of open strings and, as discussed, has an intercept an = 1/2. Physically 
this is interpreted as the exchange of quark-antiquark pairs In addition there is a Pomeron 
trajectory which has an intercept ap = ap + 1/2, and is due to the exchange of closed strings 
(interpreted as glueballs in QCD). 






Figure 8: Nucleon-parton scattering process. 



Let us now return to DIS and let us see what all this has to do with it. Let us consider 
the behaviour of the amplitude A(s,t) for large s and fixed t 

A(s,t)~ s a,t + a o. (77) 
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Let us now assume that the elastic nucleon-parton amplitude (figure 9) is described by such 
an amplitude. Kinematically, s = (P — k) 2 . If we decompose 

k 2 
k = xP- -^n + k T , (78) 

2x 

where n is a vector such that n 2 = 0, n ■ P = 1 and kx • n = kx ■ P = 0, then, with the usual 
approximations , 

s = -2k-P = ^ (79) 

x 

We see then that, at a finite value of kr, low x corresponds to the large s behaviour. On the 

other hand, for this subprocess t = 0, and we of course realize that the amplitude is directly 

related to the cross-section for parton + proton — ► anything. If the parton is a quark only the 

Reggeon will contribute. If there is mixing with partonic gluons we will have a contribution 

from both the Regge and Pomeron trajectories. After a short calculation we shall conclude 

that 

F 2 (x) ~ A P x° + A R x^ . (80) 

And, therefore that, at low values of x, 

g(x) ~ x -1 , g(x)~x~2. (81) 

These are the implications of Regge theory for the parton distribution functions. Notice 
that there is no Q 2 dependence anywhere, so the question poses itself as to which is the 
appropriate value of Q 2 to compare with. The answer is not very well defined, but it should 
correspond to the typical range of energies where Regge phenomenology is known to be valid 
in other contexts, i.e. a few GeV. In fact, a fit to the gluon PDF shows that at Q 2 = 4 GeV 2 , 
g(x) ~ x' 1 - 17 . Not bad. 

9 Low x region 

The region below 10 -2 had not been explored experimentally until very recently; a first look 
at these low-x values has been provided by the commissioning of HERA. HERA is a machine 
ideally suited for an in-depth analysis of structure functions. It should be possible to arrive 
at very low values of x (down to x ~ 10 -5 ). 

Most parametrizations have traditionally performed very poorly when extrapolated to 
the low x region. Typically they predict an increase as x — > which is lower than what is 
actually seen. The behaviour i^(^) ~ £ , with A ~ 1/2 as x — > 0, which is predicted from 
the BFKL evolution equation|£8[ seemed at some point (see e.g. [p9[| ) to stand the comparison 
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with HERA results best. However, this behaviour is still incompatible with unitarity and 
cannot hold all the way to x = either. If fact we know now that the predictions from BFKL 
cannot be trusted. This has prompted a renewed interest in trying to extract the behaviour 
at low x from conventional Altarelli-Parisi evolution. The consensus now seems to be that 
even for the low values of x analyzed at HERA there is no real evidence of any results beyond 
ordinary perturbative QCD. 

It is easy to understand why perturbative QCD must fail at some point. The expansion 
of the splitting function P(z) in powers of a s at the NLO actually resums all terms of the 
form (a s logQ 2 ) n and a n log n_1 Q 2 . Looking at (jTC| ) we see that the propagator causing the 
mass singularity is (p = £P) 

(82) 



1 
2pk 



2x 

E,Krp 



Apart from the parametric integrals, we have 



d 2 k T 

Kirp 

This is the origin of the log A 2 and, eventually, of the logQ 2 . 

The leading log™ Q 2 will thus be produced by one single region in integration 



(83) 



O Kirp 



(Jj Krp 



Ob Krp 



(84) 



with \Q\ » |fc£| > |fc£ -1 | > ...|4|. 

At sufficiently large x logarithms of 1/x necessarily appear. We have actually seen them in 
the double scaling limit. They must, at some point, spoil the predictivity of the perturbative 
expansion. One must then identify the regions in integration capable of giving rise to terms 
of the form (a s log ^) n , and eventually to a" log™ -1 \ and so on. 




Figure 9: Ordering leading to the most singular log 1/x contribution. 
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Lipatov and coworkers pq| (see alsopOfl for alternative derivations) have identified such a 
contribution. It corresponds to the diagram depicted in figure 10, more specifically to the 
region 

ki = a { P + fan + K iT , (85) 

ai » a 2 > ... > otn-i, kiT^k jT , fix < /3 2 < ... <€C A»-i- ( 86 ) 

This leads to splitting kernels similar to those of the Altarelli-Parisi equations 

F 2 (x,Q 2 ) = J d 2 k T J^C(^,Q 2 ,k T )F 2 ^,k T ), (87) 

where F 2 (£, hp) obeys the differential equation 

£ JU2& k' T ) = J d 2 k T K(k' T , k T )^F 2 (C, k T ). 



The BFKL kernel is now known to leading and subleading order. [28, 31] The leading asymp- 
totic solution is 



F(x,k T )~x- 4N[og * . (89) 



Unfortunately the corrections implied by the next-to-leading calculations are gigantic |31| . 
There is no way of doing anything useful with BFKL scaling at present. 

As previously discussed this does not seem to be a problem for HERA data since a careful 
analysis shows that — perhaps surprisingly — the data is well accounted for by ordinary 
perturbative QCD (the matter is however somewhat controversial to this date), but it will 
come the day where log 1/x corrections will be essential. The subject is thus still open. 
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